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We have performed radio-frequency dissociation spectroscopy of weakly bound ^Li2 Feshbach 
molecules using low-density samples of about 30 molecules in an optical dipole trap. Combined 
with a high magnetic field stability this allows us to resolve the discrete trap levels in the RF 
dissociation spectra. This novel technique allows the binding energy of Feshbach molecules to be 
determined with unprecedented precision. We use these measurements as an input for a fit to the 
^Li scattering potential using coupled-channel calculations. From this new potential, we determine 
the pole positions of the broad ®Li Feshbach resonances with an accuracy better than 7 x 10~* of 
the resonance widths. This eliminates the dominant uncertainty for current precision measurements 
of the equation of state of strongly interacting Fermi gases. For example, our results imply a 
corrected value for the Bertsch parameter ^ measured by Ku et al. [Science 335, 563 (2012)], which 
is ^ = 0.370(5)(8). 



PACS numbers: 

In the past few years, ultracold Fermi gases of neutral 
atoms have become important benchmark systems for 
testing theories of strongly interacting many-body sys- 
tems [1]. This success is based on two main factors. The 
first is that the physics of ultracold gases is very well ap- 
proximated by simple model Hamiltonians. These Hamil- 
tonians contain only a contact interaction, which can be 
described by a single quantity, the scattering length a. 
The second is the existence of Feshbach resonances in the 
interparticle scattering, which cause the scattering length 
to diverge to ±oo at certain magnetic field values [ ]. 
This allows tuning of the interparticle interactions by ap- 
plying a homogeneous magnetic offset field. Using such 
resonances, the properties of strongly interacting Fermi 
gases have been investigated using a number of different 
techniques, which range from radiofrequency (RF) spec- 
troscopy [ >, 4], through studies of collective oscillations 
[5, 6], to the detailed analysis of in-trap density profiles 
[7-9]. However, regardless of which technique is used, all 
such measurements depend on accurate knowledge of the 
properties of the Feshbach resonance that is used to tune 
the interactions. 

^Li atoms in the three energetically lowest Zeeman sub- 
levels of the electronic ground state (labeled |1), |2) and 
1 3) following ref. [ ]) are widely used to realize such 
strongly interacting Fermi gases. The interactions be- 
tween atoms in the three different spin states are de- 
scribed by three scattering lengths ai2, 023 and 013, 
which can all be tuned using broad Feshbach resonances 
located at magnetic fields of about 800 G [31]. The posi- 
tions of these resonances were determined by Bartenstein 
et al. with an accuracy of about IG [10], which corre- 



sponds to less than 1% of the resonance widths. However, 
recent studies of the equation of state of strongly interact- 
ing Fermi gases have reached a level of precision at which 
they are limited by the uncertainties in the resonance po- 
sitions. An important example is a measurement recently 
performed by Ku et al. [9] with the goal of measuring the 
equation of state (EoS) at the point where the scattering 
length diverges to ±cx). In this so-called unitary limit the 
scattering length drops out of the problem, leaving the 
interparticle spacing as the only remaining length scale. 
At zero temperature this has the consequence that all 
extensive quantities of the unitary Fermi gas are given 
by their values for a non-interacting system rescaled by 
a universal numerical constant ^, known as the Bertsch 
parameter. Ku et al. determined this parameter to be 
^ — 0.376 ± 0.004, providing a precision measurement 
that can serve as a test for theory. However, if the mea- 
surement is performed at a finite value of the scattering 
length, it leads to systematic errors. The error in ^ result- 
ing from the 1.5 G uncertainty in the resonance position 
determined by Bartenstein et al. is about 2% and is the 
largest error contribution [ ]. This clearly illustrates the 
necessity of a new, more accurate determination of the 
properties of the ^Li Feshbach resonances. 

In this work we determine the positions of the broad 
^Li Feshbach resonances with an accuracy of 80 mG, 
which corresponds to less than 7 x 10"''^ of the resonance 
widths. To achieve this we make use of the fact that 
every Feshbach resonance is related to a weakly bound 
dimer state. Close to the resonance the binding energy 
of the dimer is approximately related to the scattering 
length by the universal relation E-^ = fi? /ma^, where m 
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FIG. 1: RF spectra for the free-free (red dots, left axis) and 
bound-free (blue squares, right axis) transitions at a mag- 
netic field of _B = 811.139 G. The line shape of the free-free 
transition is well described by a Lorentzian (solid line). The 
bound-free spectrum is given by the wavefunction overlap be- 
tween the initial molecular state and the possible final states 
of the dissociated atoms in the trapping potential. This allows 
us to observe distinct peaks spaced by 2i/r which correspond 
to different radial trap levels. 



is the mass of one atom [ ]. Thus we can obtain infor- 
mation about the ^Li Feshbach resonances by measuring 
the binding energy of such a weakly bound dimer state 
for different values of the magnetic field. However, the 
universal relationship is not accurate enough for quan- 
titative interpretation, and in the present work we fit 
the measured binding energies to determine a new model 
interaction potential for ^Li using coupled-channel calcu- 
lations. This new potential in turn provides a{B) as a 
function of magnetic field B and allows us to characterize 
the Feshbach resonances to high precision. 

The most precise method currently available to mea- 
sure the binding energy of these dimers is RF spec- 
troscopy [10, 11]. This technique is based on applying 
an RF pulse to a gas of atoms which drives them from 
an initial hyperfine state \i) to a final state |/). For a 
sample of molecules one can either drive a transition to 
another weakly bound dimer state (bound-bound transi- 
tion) or dissociate the dimer into two free atoms (bound- 
free transition) . In either case the transition frequency is 
shifted from the free-free transition by the difference in 
the binding energies of the initial and final state. How- 
ever, the transition frequency is also affected by the dif- 
ference in the mean-field energies of the initial and final 
state. To avoid this systematic error, measurements of 
the dimer binding energy must be performed in a regime 
where the scattering length is much smaller than the in- 
terparticle spacing, i.e. na^ ^ 1. In previous experiments 
this could be achieved only for relatively small values of 
a < 2000 ao, as the experimentally achievable densities 
were limited to n ~ 10^'^ molecules/cm"^. Accordingly, 
the smallest binding energies that could be measured 
were on the order of i?b ^ 100 kHz, which resulted in 
a large uncertainty in the fitted resonance position. 



We use the techniques we have developed to prepare 
and detect few-particle systems [12] to create very dilute 
samples of molecules. This allows us to perform RF spec- 
troscopy of dimers with much smaller binding energies, 
and thus measure much closer to the resonance. We start 
from a small BEC of about 10"^ |12) molecules, trapped 
in a small-volume optical dipole trap at a magnetic field 
of 760 G. Subsequently we reduce the particle number to 
about 30 molecules by applying the spilling technique de- 
veloped in ref. [1 2]. We then superimpose a large-volume 
optical dipole trap with trap frequencies of fr = Wr/27r = 
349(3) Hz and i^ax — Wax/27r = 35(1) Hz in the radial and 
axial directions, respectively. To transfer the molecules 
into this shallow dipole trap we suddenly switch off the 
microtrap. This non-adiabatic release results in a mean 
kinetic energy per particle of > 0.4 /iK and therefore a 
final molecular density of n « 10^ molecules/cm"^. This 
corresponds to an interparticle spacing of about 10 /im. 

To measure the bound-free spectra we first perform a 
10 ms ramp from the magnetic field of 760 G at which 
we prepare the sample to the magnetic field of interest 
and wait for another 5 ms. This time is long enough for 
the magnetic offset field to stabilize, but short enough to 
avoid coUisional dissociation of a significant fraction of 
molecules. We then apply a 10 ms RF pulse to dissociate 
a fraction of the |12) molecules into free atoms in states 
|1) and 1 3). By measuring the number of these unbound 
atoms as a function of the RF frequency, we obtain spec- 
tra as shown as blue dots in figure 1. To limit saturation 
effects, we choose the pulse power such that we dissociate 
at most 30 to 40% of the molecules. From these spectra 
we extract the frequency t'bf of the bound-free transition. 

To measure the frequency of the free-free transition we 
prepare a spin-polarized Fermi gas of atoms in state |2) 
and drive the RF transition to state )3) (red dots in figure 
1). We do this before and after the molecule dissociation 
measurement and use the weighted mean i>ff of the two 
measurements (see Supplemental Material). From this 
we can also determine the magnetic field using the Breit- 
Rabi formula. 

The profile of the bound-free spectrum is determined 
by the overlap between the wavefunctions of the initial 
molecular state ^/i, and of the accessible final states ipf 
[ ] . As the RF photons carry only negligible momentum, 
there is no net momentum transfer to the system and 
therefore the RF pulse can affect only the relative motion 
of the two atoms. The transition rate between the initial 
and final state is thus proportional to \{ipi{r)\ijjf{r))f, 
where r is the distance between the two atoms. For a 
continuum of final states the resulting asymmetric line 
shape is given by the momentum distribution of the ini- 
tial molecular state [14]. In a confining potential, how- 
ever, the final states are the discrete energy levels of the 
trap and the profile is determined by the overlap between 
the molecular state and the trap states. If the experimen- 
tal resolution is insufficient to resolve the trap levels, as 
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FIG. 2: Molecular dissociation spectra at four different mag- 
netic fields. The lines show fits according to the model de- 
scribed in the text, with the solid parts indicating the range 
of the data points included in the fit. 

was the case in previous experiments [3, 4, 10, 11], the 
final state can be well described by a continuum. In our 
case the resolution is high enough to resolve the radial 
trap levels (see figure 2). 

For the initial molecular state, the long-range part of 
the wavefunction for the relative motion is well described 
by ipi{r) ~ e"*"/", where a is the scattering length. The 
final states are the levels of our optical trap, which we ap- 
proximate as harmonic oscillator levels. Since the initial 
state is symmetric, only the symmetric harmonic oscil- 
lator levels (nho=0,2,4,...) contribute. Calculating the 
wavefunction overlap results in a spectrum of delta func- 
tions of different heights located at fbf + pvr + 'Z'^ax with 



p, q positive even integer numbers. To fit our measured 
spectra, we convolute this spectral function with the line 
shape of the free-free transition, which we approximate 
by a Lorentzian with a FWHM of 122 Hz. Due to this 
finite resolution only the radial peaks are resolved. The 
free parameters of the fits are the bound-free transition 
frequency i/bf , the overall amplitude, and a small offset in 
the atom number arising from coUisional dissociation of 
molecules [l'>]. To determine i^bf we fit the lowest radial 
peak at each field (solid lines in figure 2) . For final states 
with higher energy the theoretical model for the ampli- 
tude shows only qualitative agreement with the measured 
data (dashed line). The main reasons for this are satu- 
ration effects, anharmonicity, and a small deviation from 
the assumed cylindrical symmetry of the trap. 

To deduce the molecular binding energies from the 
dissociation frequency du, we need to subtract the 
confinement-induced frequency shifts [Id] (see Supple- 
mental Material). The resulting binding energies are 
given in table L 

To fit the experimental results and extract the position 
of the broad resonance pole, we use a coupled-channel 
model similar to that of refs. [17] and [10]. The inter- 
action potentials are constructed using the short-range 
singlet potential of ref. [18] and the short-range triplet 
potential of ref. [1!)], joined at long range onto poten- 
tials based on the dispersion coefficients of ref. [20] and 
the exchange function of ref. [l!S]. The interatomic spin- 
dipolar interaction is taken to follow its long-range {r~^) 
form at all distances. The singlet and triplet scatter- 
ing lengths are adjusted by making small changes to the 
repulsive walls of the singlet and triplet potentials with 
parameters Sq and 5*1. Scattering calculations are car- 
ried out using the MOLSCAT package [21], and bound- 
state calculations using the companion package BOUND 
[22, 2--)]. MOLSCAT can converge directly on the posi- 
tions of poles and zeroes in the scattering length. We car- 
ried out least-squares fits to the new binding energy mea- 
surements described above, together with the two bound- 
bound spectroscopic frequencies of ref. [10] at 661.436 G 
and 676.090 G, the magnetic field near 527 G where the 
scattering length passes through zero [ ' ], and the po- 
sition of the narrow resonance near 543 G. The least- 
squares fits were carried out using the interactive package 
I-NoLLS [25]. 

A two-parameter fit using only Sq and 5*1 proved capa- 
ble of giving a good fit to all the experimental results ex- 
cept the position of the narrow resonance. This fit placed 
the narrow resonance about 0.12 G to high field of its ex- 
perimental position. This discrepancy could be resolved 
by introducing a third parameter in a variety of ways, 
such as scaling the exchange potential or changing the 
value of the exponent parameter /3 in the exchange po- 
tential. However, in the absence of a good theoretical 
justification for the extra parameters, and since introduc- 
ing them had little effect on the parameters of the reso- 
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magnetic field 
B (stat.) 
m Gauss 


free-free transition 
fff (stat.) 
m MMz 


bound-free transition 
fbf (stat.)(sys.) 
m MHz 


dissociation freq. 
&v (stat.)(sys.) 
m Jillz 


confinement shift 
i^cs = t'o(sys.) -f t'i(sys.) 
m kHz 


binding energy /h 
VE^, (stat.)(sys.) 
m ItHz 


811.139 (1) 
801.115 (5) 
781.057 (1) 
720.965 (1) 


81.830 115 (3) 
81.891 539 (33) 
82.019 823 (1) 
82. 452 482 (2) 


81.832 271 (7)(8) 
81.896 236 (3)(8) 
82.034 336 (6)(8) 
82.579 943 (13) (8) 


2.156 (8)(16) 
4.697 (33)(16) 
14.513 (6)(16) 
127.461 (13)(16) 


0.367(3) - 0.014(1) 
0.367(3) - 0.011(1) 
0.367(3) - 0.011(1) 
0.367(3) - 0.021(1) 


1.803 (8)(17) 
4.341 (33)(17) 
14.157 (7)(17) 
127.115 (14)(17) 



TABLE L Measured dissociation frequencies and resulting binding energies at different magnetic fields. The dissociation 
frequency &v is obtained by subtracting the free-free transition frequency from the bound- free transition frequency fbf. To 
obtain the binding energy we need to subtract the confinement induced frequency shift v^s from the dissociation frequency 5v. 
This shift has two contributions: The first is the shift vq = Vr + \vax caused by the zero-point energy of the relative motion 
of the dissociated particles. The second contribution is the interaction shift Vi caused by the confinement of the molecules and 
the interactions between the dissociated particles [16]. A complete discussion of the statistical and systematic errors can be 
found in the Supplemental Material. 



nances near 800 G, we ultimately chose a two-parameter 
fit, excluding the data point for the pole of the narrow 
resonance, as the most reliable for our purpose. To esti- 
mate the uncertainties in the pole positions and derived 
parameters, we repeated the fits using binding energies 
at the upper and lower limits of the systematic uncertain- 
ties and used the range of predictions from the various 
fits to estimate the model-dependence. 

The quality of fit and the key quantities calculated 
from the best-fit (2-parameter) potential are summarized 
in Table II. Tabulations of a{B) for the best-fit potential 
are available in [!>!]. 

With these results, the uncertainty in the positions of 
the broad ^Li Feshbach resonances is no longer a limiting 
factor for current experiments. Using our new calibration 
of a{B) it is possible to address systematic errors in re- 
cent experiments which were caused by the inaccuracy 
of the previous determination of the resonance positions. 
The most striking example of this is the determination of 
the Bertsch parameter ^ by Ku et al. [9], which was per- 
formed using a mixture of ^Li atoms in states |1) and |2) 
at a magnetic field of 834.15 G. At this field, our best- 
fit potential gives a{B) = -2.124(80) x 10^ cq and ef- 
fective range r^s = 87.03(1) Oq. The difference between 
the EoS at unitarity and the EoS measured at this finite 
value of the scattering length may be obtained by using 
Tan's contact C(a) [9, 26]. This gives a corrected value 
for the normalized zero-temperature chemical potential 
fx/Ep at unitarity, which in turn gives a revised value 
of the Bertsch parameter ^ — 0.370(5) (8) [27]. Here the 
first parenthesis denotes the statistical error, while the 
second gives the systematic uncertainty of the corrected 
value [- Vi]. 

In this work we have established a new technique to 
measure the binding energy of weakly bound molecules 
by performing trap-sideband resolved RF spectroscopy. 
By creating very dilute samples of molecules we have 
greatly reduced density-dependent shifts of the RF tran- 
sitions, which has allowed us to perform spectroscopy of 
extremely weakly bound molecules. Using these tech- 





Fit ref. [ ] 


Present fit 


Experiment 


fb,12 - ^h.V3 + m 


83 664.0(10) 


83 665.9(3) 


83 664.5(10) 


at 661.436 G 








fb,12 - I'b.lS + I'S 


83 297.3(10) 


83 297.3(3) 


83 296.6(10) 


at 676.090 G 








fb,i2 at 720.965 G 




127.115(17) 


127.115(31) 


1^,12 at 781.057 G 




14.103(26) 


14.157(24) 


1^,12 at 801.115 G 




4.342(17) 


4.341(50) 


i'b,i2 at 811.139 G 




1.828(11) 


1.803(25) 


Zero in ai2 




527.32(25) 


527.5(2) 


Narrow pole in ai2 




543.41(12) 


543.286(3) 




45.167(8) 


45.154(10) 




at 


-2140(18) 


-2113(2) 




|12) pole 


834.149 


832.18(8) 




12) A 


300 


-262.3(3) 




12) flbg 


-1405 


-1582(1) 




|13) pole 


690.43 


689.68(8) 




13) A 


122.3 


-116.6(3) 




13) flbg 


-1727 


-1770(5) 




23) pole 


811.22 


809.76(5) 




23) A 


222.3 


-200.2(5) 




23) flbg 


-1490 


-1642(5) 





TABLE II: Quality of fit between coupled-channel calcula- 
tions on the best-fit 2-parameter ^Li potential and the exper- 
iments, together with key derived quantities calculated using 
the potential. The quantities in parentheses are estimates of 
the model-dependence, including the effect of the systematic 
errors in the binding energies in Table I. All frequencies are 
given in kHz, all lengths in bohr and all magnetic fields in G. 
The A and flbg values are obtained from localHts to a{B) near 
the resonance and do not correctly reproduce the positions of 
the zeroes in a{B). 



niques we have measured the binding energy of ^Li Fesh- 
bach molecules with binding energies as low as /i x 2 kHz 
with an accuracy better than h x 50 Hz, which is a 40- 
fold improvement compared to previous measurements 
[ ] . From these binding energies we have determined 
the positions of the broad ^Li Feshbach resonances with 
an accuracy of 80 mG using a coupled-channels calcula- 
tion. This removes one of the major limiting factors for 
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precision studies of strongly interacting Fermi gases using 
ultracold atoms. 
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SUPPLEMENTAL MATERIAL 
Determination of the dissociation frequency 

To determine the dissociation frequency we measure 
the bound-free transition frequency vm and the free- 
free transition frequency vg. To check for shifts of 
the RF-transitions during the experiment we measure 
the free-free transition before (fffi) and after {i'tf2) the 
molecule dissociation measurement and use the weighted 
mean z/g of both measurements. The dissociation fre- 
quency is then given by Sv — vm — Vff. The magnetic 
field is calibrated by inserting into the Breit-Rabi- 
Formula [2S]. To determine the density dependent shift 
^I'donaity iucrcasc the particle number from 30 to about 
200 molecules. We find a shift of Sh' which is smaller 
than 50 Hz, from which we estimate a shift of less than 
0.125 Hz per particle in a linear approximation. There- 
fore the systematic uncertainty due to density effects is 
^I'dcnaity = 8 Hz for a sample of 30 molecules (60 atoms). 
All parameters involved in the determination of the dis- 
sociation frequencies are listed in TABLE HI. 



Confinement shift 

To determine the binding energy of the molecules we 
have to subtract the effects of the confining potential 
from the dissociation frequency dv. The main contribu- 
tion to the confinement shift is the zero point energy Eq 
of the relative motion of the dissociated particles in the 
trap which is given by the frequency i'q ~ I'r + \vax ~ 
367(3) Hz. In the limit of vanishing scattering length, 
a — >■ —0 , Eq is the only contribution to the shift. For 
non-zero scattering length an exact expression for the 
energy of two particles in a cigar shaped harmonic trap 
with aspect ratio 77 has been derived in ref. [Ki]. It relates 
a to the total energy i? < £'0 of the two particles 



1 

a 



1 



T{-S/2). 



(1) 



with £ — E ' 
given by 



Eq. The integral representation of J- {x) is 



dt 



Vl -e-* (1 - e-^t) 



1 



(2) 



To calculate the energy we numerically solve equation 1. 
The result for an aspect ratio of 77 = 10 is shown in fig- 
ure 3 (blue curve) . The difference to the universal bound 
state in free space (green curve, [2]) determines the con- 
finement shift. 

To obtain the confinement shift of the initial state and 
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FIG. 3: Confinement shift [29]. The green curve shows the 
energy of the universal bound state associated with a Fesh- 
bach resonance [_!]. At the point where the scattering length 
diverges at the Feshbach resonance (FR, upper panel) the uni- 
versal bound state reaches the continuum. By measuring its 
binding energy in free space (green arrow) as a function of the 
magnetic field one can determine the position of the Feshbach 
resonance. The blue curve shows the universal bound state 
in the presence of the confinement {['.'•], rj — 10). The blue 
arrow indicates an RF transition from a molecule at a > to 
two atoms with a < 0. Due to the confinement the dissocia- 
tion frequency is shifted with respect to the one in free space. 
The units are ujn — and an = 



the final state of the dissociation measurement at differ- 
ent magnetic fields we estimate the scattering length ai2 
from the corresponding dissociation frequencies. There- 
fore we first estimate the binding energy by subtracting 
uq from the dissociation frequency. Then we calculate the 
scattering length using the expression of ref. [ ] for the 
energy Egf of the bound state which considers effective 
range corrections to first order 



^2 



gf ^ 



^(ai2 - ay 



(3) 



with /i the reduced mass and with the so-called mean 
scattering length 



0.487r. 



vdw 



(4) 



where Tvdw is the range of the van-der-Waals potential. 
The scattering length at the corresponding magnetic 
fields is taken from ref. [10]. The scattering length ai2 
of the initial state and 013 of the final state of the rf- 
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magn. field 
B [GJ 


free-free 1st 
Vffl [iVlJrlz] 


free-free 2nd 
t'ff2 [MHzJ 


free-free weighted mean 
Vs [MMzJ 


bound-free transition 
Vhi [MMzJ 


dissociation freq. 

jc. n Tu 1 
Of [kHzJ 


811.139 (1) 
801.115 (5) 
781.057 (1) 
720.965 (1) 


81.830 120 (8) 
81.891 515 (3) 
82.019 822 (2) 
82. 452 484 (4) 


81.830 113 (5) 
81.891 583 (4) 
82.019 824 (3) 
82. 452 479 (5) 


81.830 115 (3) 
81.891 539 (33) 
82.019 823 (1) 
82.452 482 (2) 


81.832 271 (7)(8) 
81.896 236 (3)(8) 
82.034 336 (6)(8) 
82.579 943 (13) (8) 


2.156 (8)(16) 
4.697 (33)(16) 
14.513 (6)(16) 
127.461 (13)(16) 



TABLE III: Transition and dissociation frequencies [29] . The errors a-i.^^-^ and (Ji^jf^ of the measured free-free transition 
frequencies are the statistical errors of the fit to the free-free transition. Vfi — — is the weighted mean of um and 



U{f2 with error (t^„ — 



\ 



-('-'ffi-i'tf)^ 



The magnetic field error as is the error resulting from the statistical error a,^ 



For the bound-free transition frequency z/bf the first parenthesis gives the statistical error a^^j of a Lorentzian fit to the rising 
slope of the first peak of the spectrum. The second parenthesis gives the systematic error Ai,^^^^j of the fit, which we estimate 
by the difference between the fitted frequency using either a Lorentzian lineshape or a Gaussian lineshape to describe the 
transition peak into a single trap sideband. For the dissociation frequency Su the statistical error (first parenthesis) is obtained 
by quadratic addition of cr^jf and while the systematic error is the sum of A^^^^^^^j and the systematic uncertainty due to 
density dependent shifts A^^j^^^^j^ (see text). 



dissociation freq. 

Su [kHz] 


initial a 
ai2 [lO^bohr] 


final a 
ai3 [lO^bohr] 


initial cs shift 
Ucs-i [kHz] 


final cs shift 
i/cs-f [kHz] 


cs shift 
Ucs [kHz] 


binding energy/h 
ue^ [kHz] 


2.156 (8)(16) 
4.697 (33)(16) 
14.513 (6)(16) 
127.461 (13)(16) 


18.34 
11.80 
6.54 
2.20 


-3.54(1) 
-3.69(2) 
-4.10(3) 
-8.71(22) 


0.006 
0.002 
0.002 
0.000 


0.359(1) 
0.358(1) 
0.357(1) 
0.346(1) 


0.353 (3)(1) 
0.356 (3)(1) 
0.356 (3)(1) 
0.346 (3)(1) 


1.803 (8)(17) (25) 
4.341 (33)(17) (50) 
14.157 (7)(17) (24) 
127.115 (14)(17) (31) 



TABLE IV: Dissociation frequencies and binding energies [29]. The scattering length ai2 is calculated using equation 3 
and ai3 is determined from ai3(_B) of ref. [''.']. Its systematic error results from the 1 G uncertainty of the |13) pole in ref [ ]. 
The confinement shift of the initial and final state of the rf-transition is calculated using equation 1. The error is the propagated 
error of aia. The difference of both shifts determines the total confinement shift i/cs- The first parenthesis states the statistical 
error of the zero point energy of the relative motion given by the SEM of the radial trap frequency which is determined from 
the separation of the sideband peaks in the dissociation spectra. The second parenthesis gives the systematic error of the 
confinement shift of the final state. The quadratic addition of both errors determines A^^^,. The binding energy Eb is calculated 
from the difference between the measured dissociation frequency and the confinement shift. The first parenthesis gives the 
statistical error cji/^^ (see TABLE III). The second parenthesis gives the systemtic error A^/^^ = A,^^^^^j+A,^^^^^.^^+Ai,^^. The 
third parenthesis gives the sum of the statistical and systematic error. 



dissociation measurement and the corresponding confine- 
ment shifts are Hsted in TABLE IV. By subtracting the 
total confinement shift from the dissociation frequency 
we obtain the binding energy of a molecule in free space. 



